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Heisenberg-Weyl algebra is called a TAA algebra after Tekin, Aydin, 
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to use a matrix approach to implement analytic representations of the 
Heisenberg-Weyl algebra in univariate and multivariate settings. We 
provide examples for the univariate case. Krawtchouk polynomials are 
presented in detail, including a review of Krawtchouk polynomials that 
illustrates some curious properties of the Heisenberg-Weyl algebra, as 
well as presenting an approach to computing Krawtchouk expansions. 
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how to implement in finite terms Rota's "finite operator calculus" . 
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1 Introduction 



Ever since Weyl's "Group Theory and Quantum Mechanics", the question 
of finite-dimensional representations of the canonical commutation relations 
has been of interest. The simplest formulation, two operators A,B on a, 
Hilbert space obeying 

[A, B] = cl 

for a non-zero scalar c, 1 being the identity, is quickly dispatched on a finite- 
dimensional space 7i by taking traces, the left side yielding zero, the right 
giving c dimTi. 

We illustrate with the interesting phenomenon of "Krawtchouk ghosts" , where 
the Krawtchouk polynomials, polynomials orthogonal with respect to a bi- 
nomial distribution, are an infinite family of polynomials that carry a rep- 
resentation of the HW commutation relations, but span a finite-dimensional 
Hilbert space, all but a finite number of the polynomials having zero norm. It 
turns out that the Krawtchouk polynomials in each dimension are a basis for 
corresponding representations of su(2), used ubiquitously in current work on 
quantum information. However, they do not appear explicitly yet in current 
work on quantum information. 

Mathematical physicists have used Krawtchouk polynomials to develop finite- 
dimensional quantum mechanics. See Lorente [8, 9] for work related to this 
approach. Atakishiyev & Wolf [1, 2] have used Krawtchouk polynomials to 
formulate finite-dimensional wave functions to replace Gaussian packets, e.g. 
in optics. Santhanam [13] has discussed the difficulties as well as possibilities 
of finite-dimensional quantum mechanics. In other domains, Krawtchouk 
expansions have been used in image analysis showing results superior to 
many commonly used techniques [17]. 

The approach here will show how up to a given order one may produce 
representations of the HW algebra on spaces of polynomials that we call 
"analytic representations" , whereby corresponding to any function analytic 
in a neighborhood of the origin in C, one produces a sequence of canonical 
polynomials, that are closely related to many of the polynomial sequences 
arising in physics, probability theory, combinatorics, etc. [5, 11, 12]. 

Another approach to finite-dimensional representations of the canonical com- 
mutation relations is via replacing the identity by an operator E with the 
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defining relations 

[A,B]=E, [A,E] = [B,E] = Q (1) 

whereby E is central. This is ilhistrativc of the inapplicability of Schur's 
Lemma where a central element must be a multiple of identity as for finite- 
dimensional representations of unitary groups. The finite-dimensional rep- 
resentations of the relations (1) have been found using a diagrammatic ap- 
proach in [7]. 

There are three main aspects of this work: 

1. We make the observation that the TAA algebra is the appropriate for- 
mulation for abstracting the structure of the truncated Hciscnbcrg- 
Weyl algebra. We mention related work in the area of discrete/finite 
quantum mechanics. 

2. We show that for analytic representations of the HW algebra, compu- 
tations using truncated HW can be done order-by-order numerically, 
avoiding the necessity for symbolic computations. 

3. We show some properties of Krawtchouk polynomials that are of inter- 
est with regards to HW representations. 

The algebra we call here "TAA algebra" was presented by Tekin, Ay din, and 
Arik in [14]. It is generated by an operator a and its adjoint a*. Setting 
V = a*a, the commutation rule defining the algebra is [a, u] = a. Taking 
adjoints gives the complementary rule [u, a*] = a*. In other words, we assume 
that u and a generate the two-dimensional Lie algebra of the affine group. 
Writing it out 

aa*a — a*aa — a (2) 

we see that the commutation rule aa* — a*a — 1 oi the Heisenberg-Weyl 
algebra has been multiplied by a on the right. This modification is enough to 
yield finite-dimensional representations, including the truncated Heisenberg- 
Wcyl algebra (HW-algebra) given by the operators X =multiplication by x 
and D = d/dx acting on polynomials of a given bounded degree. Writing 
matrices for these operators we will see that they obey equation (2). 
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This article may be thought of as reahzing Rota's idea of "Finite Operator 
Calculus" [11, 12] in a truly finite way. There is work of P. R. Vein along 
similar lines [15, 16]. The main feature here is that the operator calculus 
is done on finite- dimensional spaces and can be carried out explicitly using 
matrices. The approach in this paper is based on algebraic properties of 
the operators and includes a formulation for the multivariable case. The 
one- variable case is dual to that presented in [4] . 



2 Orthofermion formulation 

Here we recall the orthofermion approach of [14] . Start with a set of operators 
{ci, . . . , Cp}, with p a positive integer and form the star-algebra generated by 
the {cj} modulo the following relations 

CiCj — 

p 

CiC* + 6ij J2 c*kCk = Sij 1 (3) 
fe=i 

p 

where 1 is the identity operator. Setting U = 1 — J2 c*kCk (as in [10]) we can 

fe=i 

write this last relation as 

CiC* = 6ij n 

where one readily shows from the defining relations (3) that 11^ = 11, i.e., 11 
is a projection as suggested by the notation. 

It follows from the defining relations, that Uck — Ck and from the second 
relation of eq. (3) we have the useful relation 

CiC*Ck = 5ij Ck (4) 

Within the orthofermion algebra, following [14], modifying slightly their for- 
mulation, we set 

p 

k=2 
P 

«^ = c*i + Y.clck-i 

k=2 
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Using equation (4), we get 

aa' — a' a = 1 — (p + 1) c*Cp 

which then yields the relation corresponding to equation (2) of the TAA 
algebra. 



3 Calculus with matrices 

Notation. We will denote the matrix corresponding to an operator by using 
a "symbol. Thus, X is the matrix corresponding to the operator X, etc. 

Restricting the differentiation operator to the finite-dimensional space of 
polynomials of degree less than or equal to p is no problem. Use the standard 
basis {1, x,x'^, . . . , x^}. For p = 4, we have 





f 
















2 

















3 

















4 


Vo 











0/ 



with the extension to general p following the same pattern. However, mul- 
tiplication by X must be cut off. If we define X x^ — x^'^^ for i < p and 
Xx^' = 0, we no longer have the relation DX — XD — 1. Instead, we have 
the TAA relation 

DXD - XDD = D 
The matrix of X has the form, for p = 4, 



/O 0\ 

1 

X^ 10 

10 

Vo 1 0/ 



Note that X^^^ = 0. To keep in line with the powers of x, we label the 
basis elements starting from 0. So let denote the column vector with the 
only nonzero entry equal to 1 in the {k + 1)^* position. The vacuum state is 
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Q — eo, satisfying DQ — 0. And X'^Q — e^, ior 1 < k < p. As expected, 
these are raising and lowering operators satisfying 

= e^+i 6kp , -Defe = k e^.i 

where 9ij — 1 if i < j, zero otherwise. 

With the inner product (e„, e^) = Snmnl , we indeed have D* — X. 

Let Eij denote the standard unit matrices with all but one entry equal to 
zero, {Eij)ki = SikSji, 1 < k, I < p + 1. The connection with orthofermions 
is given by the x matrix realization 

for 1 < i < p. The orthofermion relations hold and particularly for this 
realization 

Note that n = En and that the star-algebra generated by the Cj is the full 
matrix algebra. 

As long as X never multiplies the power x^, the matrix implementation agrees 
with usual calculus. The TAA relation formulates this algebraically. 

The following theorem shows that D and X not only do not generate a 
Heisenberg algebra, but, in fact, are as far as possible from doing so. 

Theorem 3.1 Forp > 0, let D and X be (p+l)x (p+1) matrices defined by 
p ^ p 

D — ^^k Ekk+i, = X! ^k+ik- Then the Lie algebra generated by {X, D} 

k=l k=l 

is sl{p + 1) . 

Proof: For convenience set n = p + 1. First we have 

[A^] =-pKn + E^fefe 
k=l 
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Set ^1 = X, 771 = b, and Hi^ H. For 2 < A; < n, let = //(adX)^ and 
77jfc = {eidb)''H, where (adA)5 = [A, S] and A{adB) = Then it is 

easily checked by induction that 

= -nEnn-k+i and Vk = akd 

for nonzero constants afc, the f denoting matrix transpose. Thus, we obtain 
Ein and Eni for 1 < i < p. Noting that [E^n, Enj] = E^j if i 7^ j, we have all 
of the off-diagonal £"s. And 

Hk — [Vk, Cfe] — —na^ {En-k+ln-k+l — Enn) 

fill out the Cartan elements of sl{n). □ 
3.1 Examples 

Here we look at some important operators. 

Example 3.2 The number operator is XD. For p = 4, we have 



( ° 











0\ 





1 

















2 

















3 





u 











4 j 



This operator multiplies by n, for < n < 4. In general, we have 

p 

XD —'^n En+l n+l 

n=0 

which multiplies by n, for < n < p. 

Example 3.3 The Hermite polynomials, occurring in oscillator wave func- 
tions, are eigenfunctions of the Ornstein-Uhlenbeck operator, XD — tD^, 
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t > 0, which for p = 4 takes the form 
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The eigenvector for each eigenvalue A = 0, 1, 2, 3, 4 gives the coefficients of the 
corresponding polynomial Hx{x, t). The family of polynomials {Hx{x, t)}AeN 
provide an orthogonal basis for with respect to the Gaussian measure with 
mean zero and variance t. 

Example 3.4 The translation operatorTt — e*-^ acts on functions as e*^f{x) — 
f{x + t). Forp = 4, 

f 1 t t^ t^ t^\ 
1 2t 3t^ At^ 
f,= 1 3t 6t^ 

1 4t 

V 1 y 

generally, with cohimns given by binomial coefficients times powers of t, 
corresponding to the action x ^ x + t on the basis polynomials x^ . The 
matrix can be computed as the exponential of tD defined as a power 
series: l + tD + t'^D'^/2\ H . 

Example 3.5 The Gegenbauer polynomials satisfy 

[(XD + a)^ - D^]C^(x) = (n + a)' C^{x) 

see, e.g. [3]. Thus we have the Gegenbauer operator, Ga — {XD + a)^ — D'^, 
which for p — A takes the form 
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where the spectrum is evident along the diagonaL Up to order p, one obtains 
the Gegenbauer polynomials with coefficients given by the eigenvectors of 

3.2 Multivariable calculus with matrices 

Here we extend to N variables. For matrices, A, B, the tensor product A^B 
denotes the Kronecker product of the two matrices. That is, if A is n x n, 
and B is m X m, then A ^ B is nm x nm with entries formed by replacing 
each entry in A with the block matrix aijB. For products of more than 
two matrices, we conventionally associate to the left. 

For a fixed p, we have (p + 1) x (p + 1) matrices 1) and X. Let / denote the 
(p + 1) X (p + 1) identity matrix. Then we set 

Dj = I(g)I(g)---(g)D(g)I---(g)I (D inthe spot) 

Xj = / (g) / • • • (g) X (g) / • • • (g) / (X in the j*"" spot) 

Then Dj and Xj will satisfy the TAA relations while [Dj,Xi] — [Xj,Xi\ — 

4 Analytic representations of the HW-algebra. 
Canonical polynomials 

Now we would like to discuss analytic representations of the HW-algcbra. 
These are infinite-dimensional representations in the sense that they act 
on a basis for the vector space of polynomials in a given set of variables 
{xi,X2,---,xn}. Basic to our approach is the use of canonical variables 
which are functions of X and D obeying the HW relations on an infinite- 
dimensional space, which restricts to the TAA relation on spaces of polyno- 
mials in a; of a given bounded degree. 

Let us review the basic construction and notations for the general, multivari- 
able, case. 

Notation. We use the convention of summing over repeated Greek indices, 
irrespective of position. 
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Given V: ^ C^, V{z) = {Vi{zi, ...,zn),..., VV(-Zi, ■ ■ ■ , zn)) holomor- 
phic in a neighborhood of the origin, satisfying V^(0) = 0, we construct an 
associated abehan family of dual vector fields. Corresponding to the oper- 
ators Xi of multiplication by Xi, we have the partial differentiation oper- 
ators, Di. In this context, a function of x — {xi, . . . ,xn), f{x), is iden- 
tified with the operator of multiphcation by f{X) acting on the 
vacuum state 1, with Dil = 0, for all 1 < i < A^. We define operators 
V{D) = (Vi^Di, . . . , Dn), . . . ,Vn{Di, . . . , Dn)). These are our canonical 
lowering operators, corresponding to differentiation. 

Denoting the Jacobian j by V'{z), let W{z) = {V'{z))-\ be the mverse 

(matrix inverse) Jacobian. Then the boson commutation relations extend to 
dVi 

[Vi{D),Xj] — -wj^- Now define the operators 



These are our canonical raising operators, corresponding to multiplication 
by Xi. We have 

Thus, the canonical system of raising and lowering operators is {K,}, {V^}, 
1 < < The essential feature, which has to be checked, is that, 
[Fj, Yj\ = \Vi, Vj\ = 0. Notice that exchanging D with X is a formal Fourier 

d 

transformation and turns the variables into the vector fields Yi = W{x)^i — 
Thus, the 1^ are dual vector Gelds [5] . 



dx^ 



Notation. We complement the standard notations used along with V and 
W, letting U denote the inverse function to V. I.e., U o V — V o U — id. 
Explicitly: U{V{z)) = z. 

Observe that since W — V'~^, we have W{z) — U'(y{z)). In other words, 
converting from zioV acting on functions of the canonical variables Y^, we 
have the recurrence relation 

X^YU'{V)-^ . 
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Using multi-index notation, n = (ni, . . . , tin), v'^ — v^'^V2^ ■ ■ ■ , the main 
formula (cf. [5, p. 185, eq. (1)]) is 



exp(v^F^) 1 = expXu,Uf,{v) = X) -r Vnix) 

n>0 

This expansion defines the canonical polynomials: yn{x) = ¥"■!. 
4.1 Canonical Appell systems 

An Appcll system, in one variable is a system of polynomials provid- 

ing a basis for the vector space of polynomials with deg /i„ = n, n = 0, 1, 2 . . ., 
such that Dhn = nhn-i- Defining the raising operator R by Rhn = hn+i, 
we have [D,R] — 1, thus a representation of the HW-algebra. (See [6, v. 3, 
Ch. 1] for further elaborations.) 

Introducing a Hamiltonian H{z), in this context the only requirement being 
analyticity in a neighborhood of the origin in C^, we have the time-evolution 



exp {-tH{D)) e^^(") = ^xU(v)-tHmv}) = ^ ^ y^{x, t) (5) 

n>0 ^• 

An Appell system of polynomials has a generating function of the form 
exp {xz - tH{z)) = X] n" ^) 



n>0 

For the canonical Appell system we have 



exp {xz - tH{z)) = J2 Vnix, t) (6) 

and we recover (5) via the inversion z = U{v), which we interpret as changing 
to canonical variables. 

Observe that each of the polynomials yn{x,t) is a solution of the evolution 
equation 
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5 Canonical calculus with matrices 



First consider the case iV = 1. We have a function V{z) analytic in a 
neighborhood of the origin in C, normahzed to V{0) — 0, V^'(O) 7^ 0. Let 
W{z) — 1/V'{z) have the Taylor expansion 

W{z) =wq + wiz H h Wkz'' H 

The corresponding canonical variable is y = XW{D), satisfying y] = 

1. The canonical basis polynomials are yn{x) — V^l, n > 0. 

Fix the order p. Let W — W{D). Then we employ the algebra generated 
by the operators V = V{D) and Y = XW. Note, e.g., that since pP+^ = 0, 
the operators V and W are polynomials in D. Similarly, since X^"*"^ = 0, 
the polynomials yn{X) are truncated if n > p. However, for n < p, the 
correspondence between the polynomials yn{x) and vectors y„ = |/„(X)eo is 
exact. Namely, the vector y„ gives the coefficients of the polynomial |/„(a;). 
The reason this works is that up to order p, the operator X never acts on a 
power of X greater than p. 

5.1 Examples 

Example 5.1 A simple example to illustrate the construction is given by 

V{z) = e'-1, U{v) = log{l + v) 

so W{z) = e-', Y = Xe-^. The relation X = YU'{V)-^ reads X = Y + YV 
or xyn = yn+i + nyn yielding the recurrence 

yn+i = (x - n)yn 

for n > 0. From yo = 1, we easily calculate 

yn{x) = x{x - 1) ■ ■ -{x - n + 1) . 
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For p — 4:, with Y 
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we get 
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with the first column giving the coefficients of the corresponding polynomial 
Un- where, since the leading coefficient equals one, we can see the truncation 
beginning in this last. 



Example 5.2 Another interesting example is the Gaussian with drift a > 0, 



V{z) = az- , U{v) = a- Va^ _ 2v 
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the minus sign taken in U{v) to have C/(0) = 0. Then W{z) 



a — z 



-, and 
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Powers of Y yield the canonical polynomials, the first few of which are 
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y4 = 15 — + 15 — + 6 — + — 
a' a° a'^ a* 

2 3 4 

/Tf» ry^ ry'-' ry 

1/5 = 105 — + 105 — + 45 — + lo- 
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These are a scaled variation of Bessel polynomials. 

In this case U'{V) ^ — ail ) . Thus, expanding and rearrangin] 

the relation X = YU'(V)-\ 

Y-X y(— -— -— -— -— — — — — 

~ 1^^^2^^2^^8^^8^'^T6^^16^^" 

which translates to 



n nin — 1) n(n — l)(n — 2) 

a yn+l = Xyn+ -Vn^ ^ Vn-l H 7^ yn-2 + ■■■ 

n 



2q;5 



" n\ (2A';-3)!! 
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Example 5.3 Our final example in this section, involves the LambcrtW 
function, which we denote W to avoid confusion with our W. Take V{z) = 
ze-^ [6, V. 1, p. 110]. Then U{v) = -W{-v). We find Y = Xe^{I - D)-^ 
and with p = 7 the corresponding matrix 



Y = 
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One can show that y„ = x{x + n)" ^ and that the relation X — YU'{V) 
leads to the recurrence 

yn+i = {x + 2n)yn + X) L _^ ^ 1 fc'^' Vn-k ■ 



5.2 Matrix expansions 

Using the unit matrices Eij of size (p + 1) x (p + 1), we can write formulas 
for the basic operators. We have used in Theorem 3.1 the expressions 

p ^ p 

D = Y^kEkk+i and X = Y^Ek+ik 

k=l k=l 

Then induction yields 

l<k<p 
l<k+j<p+l 

with {k)j — k{k + 1) • • • (A; + j — 1) denoting the rising factorial. Multiplying 
by X gives 

y= ^ Ek+ik+j{k)jWj 

l<k<p 
l<k+j<p+l 
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For > 1, the matrices for Dj and Xj provide the operators Yj — X^W^i(£)) 
as matrices. Repeated multipUcation on the vacuum vector Bq yields exactly 
the coefficients of the polynomials y„ up to order p, that is, no variable Xi 
appears to a power higher than p. 

6 Krawtchouk polynomials 

Note. In this section N will indicate a discrete time parameter taking values 
N — 0,1,2, We will be working with polynomials in one variable x. 

The Krawtchouk polynomials occur as polynomials orthogonal with respect 
to the binomial distribution. Here we take the distribution of the sum of N 
independent Bernoulli random variables taking values ±1 each with proba- 
bility 1/2. We start with the generating function 

G{v; X, N) = {1 + (1 - = ^ ^ K^{x, N) (7) 

where for x the position of the random walk after N steps, (A^ + x)/2 is the 
number of positive jumps and {N — x)/2 the number of negative jumps. To 
see this in the form of a canonical Appell system, write 

G{v;x,N) = e^^^'')-*^^^^^)) = {1^^ {1 - vY^' 
where we identify t ■(-^ N, 

U{v) = i log and H{U{v)) = ~ log(l - v^) 



so that 

V{z) — tanh z and H{z) — log cosh z 
with W{z) — 1/V'{z) — cosh^ z. Now verify that 




G{x,v;N) = e-^^(^)G(a;,^;;0) = (sech D)^e^^(^) 
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appropriately for the canonical Appell system. 

The recurrence relation for {i^„} is derived as follows. The recurrence rela- 
tion for the canonical polynomials with N — Ois X — YU'{V)~^ — Y{l—V^). 
Let yn{x) — Kn{x, 0). Then Kn{x, N) — (sech D)^ yn{x) . The recurrence for 
Vn is 

xyn = Vn+i - n{n - 1) (8) 

Now use the relation [f{D),x] ~ f'{D) for a function / analytic in a neigh- 
borhood of the origin to get 

[(sechD)^,x] = -A^(sechD)^-^(sechD tanhD) 

= -Ar(sechL>)^tanhL> = -A^e-^^(^V(L>) (9) 

thus, applying (sechD)^ to (8) yields 

xKn- NnKn_i = Kn+i - n{n - l)Kn-i 

or 

xKn^ Kn+1 + n{N-n+ (10) 
with K_i = 0, Ko = 1. 

6.1 Krawtchouk ghosts 

The raising operator for the polynomials when A'" = 0, is 

Y ^XW(D)^xcosh'^D 

Since Kn = (sechD)-^7/„, we have the raising operator R, RK^ = K^+i, 
given by 

R = (sechD)^r(sechD)-^ = (sechD)^F(coshD)^ 

= (sech D)^x(cosh = x cosh^ D-N sinh D cosh D 

using (9). Acting on polynomials, with V{D) — tanhD, we have the com- 
mutation relation 

[V{D),R]^1 

So we have a representation of the HW-algebra with the commutator equal 
to the identity. 
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Now notice that ifwesetn — N + 1 and n — N + 2in (10), we get 

xKn+1 = Kn+2 and xKn+2 = - {N + 2)Kn+i 

so that for n > N, all have i^^Ar+i as a common factor. Expanding the 
binomials in (7), we have 




The random walk can land on any of the points N, N — 2, . . . ,2 — N, —N in 
steps. Let x = N — 2j , where < j < N, j being the number of negative 
jumps of the random walk. Substituting into the generating function (7): 

G{v; N - 23, TV) = (1 + v)''-^ (1 - vY = E ^ ^n(A^ - 2i, N) 

The left side is a polynomial in v of degree A^, so that all of the coefficients 
beyond N vanish identically. So, Kn{N - 2j, N) = for < j < N, n > N. 
That is, on the support of the binomial distribution, i^Tjv+i vanishes. Thus 

= ll-fCAr+ill^ = ||iriv+2||^ — . . . — IliTjv+fell^ = • • • 

for k > 1, i.e., the L^-norms of all of the polynomials vanish for n > N. 
So in the Hilbert space, these are zero. Prom the HW-algebra point of view, 
these are thus "ghost states" . 

On the other hand, if we consider the lowering operator L satisfying 

LKr^ = n{N - n + l)-fr„_i 

We have the commutation relations 

[L, i?] = A , [R, A] = 2i? , [A, L] = 2L 

which give a representation of sl(2). In fact, for > 0, we recover the 
irreducible representations of su(2). 
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6.2 Krawtchouk calculus with matrices 

In the Krawtchouk Hilbert space, we consider only functions on the spectrum 
of the operator X, namely the finite set of points N, N — 2, . . . ,2 — N, —N. 

Wc have an (iV + l)-dimcnsional space spanned by the polynomials Kn{x, N), 
< n < N of bounded degree. So in this space, the matrix representations 
X and D of order (iV + 1) x (iV + 1) will give exact results. 

Fix an order p and time parameter N. Take p — N to get the full basis 
for time N. First, construct Y = Xcosh D, using e"^^. Compute F", for 
< n < p. The first column of are the coefficients of the "time-zero" 
polynomial y„. Let Sn = (sechi))^, using sechD — (cosh£))~^. Then, up 
to order p, for any N >0, 

kn{N) = SnY„ 
has first column the coefficients of Kn{x, N). 



Example 6.1 We illustrate with N ^ 5. We find 



cosh D 



and 





/ 1 





1 





1 














1 
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I 
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) 




/ 1 





-1 
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-3 







25 
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/ 
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The raising operator at time-zero is 



( ° 



















1 
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1 
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40 
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20 




I 











1 





J 



Taking p = 5, we collect 
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/ 
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496 
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120 
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168 





2720 










Y' = 










































1 





18 





280 






1 







20 





504 


















1 





32 












1 







38 





680 


















1 





20 


J 




V 







1 





32 





) 










/ 















\ 






/ 




























24 





640 





10880 








24 





688 







18752 











8 


184 




4800 
















184 





5904 







103360 




Y' = 












































20 





624 





10880 








20 





664 







18528 
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40 




1144 
















40 





1384 







24480 






i,0 1 





38 





680 / 
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1 





40 







1144 





/ 



corresponding to the time-zero polynomials |/n, < n < 5. These work for 
any N . For example, for K^{x, 3), compute 



9 





768 





23352 








9 





1294 





25160 


-10 





-536 





-15792 






-10 -516 -9520 

40 1144 

1 38 680 y 

20 



(sechL>)''y^ 




1 

V 



{sech D)^Y^ 






2480 





88120 





149 





7728 





227032 





-1016 





-35616 





-30 





-1336 





-38672 





40 





1384 





1 





40 





1144 



or K^{x, 3) = - lOx^ + 9 = - l^){x^ - 3^), vanishing at x = ±1, ±3. 

And K4{x,5) corresponds to 

1564000 \ 


-631040 


24480 
J 

Note that, in fact, one only needs the first cohimn of for each n, then 
applying appropriate powers of sech D produces the polynomials for > 0. 

Here are the polynomials for A" = 5: 

Kq = 1, Ki = X, 

K2^x^- 5, Ks^x^- 13x, 

Ki^x^- 22^2 + 45, K^^x^- 30^^ + 149x 

and 

Kq^x^ - 2>bx^ + 259a;^ - 225 = {x^ - l^){x'^ - 2,'^){x^ - 5^) 
accordingly. 



6.3 Krawtchouk expansions 

The matrix method can be used to compute Krawtchouk expansions for 
functions / defined on the spectrum -A^, 2 - A^, . . . , A^ - 2, A^ for A^ > 0. 
Start with 

G{v; X, N) = ( ) (1 - vy/' = ^ ^ K^{x, N) 

Substitute v = V{z) = tanh^; and rearrange to get, cf. (6), 

= (cosh ^^^^ A^) 



n>0 
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Replacing z by Dg — d/ds, apply both sides to a function f{s): 

e^^V(5) = f{s + x) = (cosh Dgf Y: ^^^^^^ K^{x, 7V)/(s) 



n>0 



ni 



Letting s = 0, thinking of / as a function of x instead of s, we can replace 
Ds by our usual D — d/dx, to get 

Kn{x,N) 



E 

0<n<7V 



(coshL»)^(tanhL»)V(0) 



In other words, the coefficients of the Krawtchouk expansion of f{x) are 
given by 

f{n) = ^ (coshL»)^(tanhL»)"/(0) = — (coshL>)^-"(sinhL»)V(0) 

If f{x) is a polynomial in x, we can construct a vector from its coeffi- 
cients, apply the matrices (cosh£))-^(tanh£))", rescahng by n\, and find its 
Krawtchouk expansion. 

Example 6.2 Here are the matrices (cosh£>)^(tanhi))"/n! for N — 5, 
0<n<5. Forn = 0,1,2: 



/ 1 5 65 \ 

1 15 325 

1 30 

1 50 

1 

V 1 J 
For n = 3, 4, 5: 

/ 1 30 \ 

4 

10 





Vooooo oy 



/ 1 13 241 \ 

2 52 

3 130 

4 

5 

VOOOOO oy 



/ 1 22 \ 

3 110 

6 







V 





10 


y 



/ 1 \ 

5 







Vooooooy 



/ 1 \ 









vooooooy 
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Let f{x) — x'^ + 2x^ — x'^ + 5x. One checks by hand, with N — 5, that 



f = K4 + 2K3 + 2IK2 + 31Ki + 60 

Form the column vector, transpose of [0, 5, —1, 2, 1, 0] and apply each of the 
above matrices. Stacking these as row vectors we have the matrix 

/ 60 35 29 2 1 \ 

31 50 6 4 

21 6 6 

2 4 

1 

\ 000/ 



whose first column gives the Krawtchouk coefficients. 

Alternatively, start with the top row of (coshZ))^. 
recursively, for 1 < n < A", 

Un = Un-i tanh D/n 



Call it Uq. Generate 



and construct the matrix Y with rows yo, . . . ,yn, ■ ■ ■ ,yN- Applying to the 
vector of coefficients of a polynomial gives the coefficients of the Krawtchouk 
expansion of that polynomial. In particular, applying to the vector of coeffi- 
cients of a Krawtchouk polynomial gives a standard basis vector, zero except 
for a single entry equal to one. That is, is a matrix whose columns 
are the coefficients of the corresponding Krawtchouk polynomials. For the 
example above, we have 



Y 



/ 1 





5 





65 










1 





13 





241 
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22 
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30 
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V 
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J 
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and 



/ 1 





-5 





45 





\ 





1 





-13 





149 










1 





-22 
















1 





-30 
















1 







VO 














1 


/ 



7 Summary and prospects 

In this article, the one-variable case of the matrix approach has been pre- 
sented in some detail, along with the basic theory for the multivariate case. 
The TAA algebra conveniently replaces the HW algebra in the finite-dimensional 
setting. The connection with orthofermions is interesting and clarifies the 
underlying structure. 

The approach here allows for doing calculations in finite, exact terms. There 
are many possible applications including image analysis, optics, as well as 
quantum information. 
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